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including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

08  DEC  2005 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2005  to  00-00-2005 


4.  TITLE  AND  SUBTITLE 

HY COM  Code  Development 


6.  AUTHOR(S) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Stennis  Space  Center, MS, 39529 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 

9th  HYCOM  Consortium  Meeting,  Dec  6-8, 2005,  RSMAS,  Miami,  FL 


14.  ABSTRACT 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OE  PAGES 

Same  as 

18 

Report  (SAR) 

19a.  NAME  OE 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


HYCOM  2.2  (I) 


Maintain  all  features  of  HYCOM  2.1 

o  Orthogonal  curvilinear  grids 

o  Can  emulate  Z  or  Sigma  or  Sigma-Z  models 

o  It  is  “Arbitrary  Lagrangian-Eulerian”,  see: 
Adcroft  and  Hallberg,  O.  Modelling  1 1  224-233. 

o  Explicit  support  for  1-D  and  2-D  domains 

o  KPP  or  Kraus-Turner  or 
Mellor-Yamada  2.5  or  Price-Weller-Pinkel 

o  Rivers  as  bogused  surface  precipitation 

o  Multiple  tracers 

o  Off-line  one-way  nesting 

o  Scalability  via  OpenMP  or  MPI  or  both 

o  Bit-for-bit  multi-cpu  reproducibility 


HYCOM  2.2  (Ila) 

Alternative  scalar  advection  techniques 

o  Donor  Cell,  FCT  (2nd  and  4th  order),  MPDATA 
o  FCT2  replaces  MPDATA  as  standard  scheme 

Vertical  coordinate  changes 

o  Vertical  remapping  uses  PLM  for  fixed 
coordinate  layers 

o  Thin  deep  iso-pycnal  layers 

o  Stability  from  locally  referenced  potential  density 

o  Spatially  varying  layer  target  densities 

o  Different  isopycnal  layers  in  semi-enclosed  seas 


HYCOM  2.2  (Ilb) 

Special  halo  exchange  for  tripole  global  grid 

o  Arctic  dipole  patch  on  standard  Mercator  globe 
o  Logically  rectangular  domain 
o  Two  halves  of  top  edge  “fold”  together 
o  V-velocity  changes  sign  across  the  fold 

Improved  thermobaricity 

o  No  single  reference  state  is  approriate 
for  the  global  ocean 

o  Hallberg,  Ocean  Modelling,  8,  279-300 

o  Use  a  linear  combination  of  two  out  of  three 
reference  states 

o  Atlantic  (3°C,  35.0  psu) 

o  Arctic/Antarctic  (0°C,  34.5  psu) 
o  Mediterranean  (13°C,  38.5  psu) 
o  Most  locations  use  just  one  reference  state 

o  Linear  combinations  allow  smooth  transition 
between  states 

•  Do  this  in  shallow  water  if  possible 


1/12°  GLOBAL  THERMOBARIC  REFERENCE 

STATE 


TBARIC  MAP  (1=Arctic,2=Atlantic,3=Med.) 


HYCOM  2.2  (Ilia) 

Mixed  layer  changes 

o  GISS  mixed  layer  model 
o  KPP  bottom  boundary  layer 
o  KPP  tuning 

Atmospheric  forcing  changes 

o  Option  to  input  ustar  fields 
o  Otherwise  calculated  from  wind  stress  or  speed 
o  Can  relax  to  observed  SST  fields 
o  Improved  COARE  3.0  bulk  exchange  coefficients 
o  Black-body  correction  to  longwave  flux 
o  Climatological  heat  flux  offset, 

Q  =  iQsw  ~  Qiw)  +  (Qi  +  Qs)  +  Qc 

o  Qc  is  constant  in  time 

•  Typically  based  on  the  model’s 
climatological  SST  error,  times 
(say)  -45  Wm-^QC 


deg  C 


ADDED  AFTER  FIVE  YEARS 
GLOBAL  MEAN  SST  and  SSS 


GLBaO.72,  Basin-Wide  Mean  SST 


GLBaO.72,  Basin-Wide  Mean  SSS 


LULU 


ADDED  AFTER  FIVE  YEARS 
GLOBAL  MEAN  SST  and  SSH 


GLBaO.72,  Basin-Wide  Mean  SST 


GLBaO.72,  Basin-Wide  Mean  SSH 


HYCOM  2.2  (Illb) 

Improved  support  for  rivers 

o  Still  bogused  surface  precipitation 
o  High  frequency  inter-annual  river  flow  allowed 
o  Add  it  to  atmospheric  precip,  off-line 
o  Instead  of  monthly  climatology, 
or  in-addition  to  it  (flow  anomalies) 

o  Better  control  of  low  salinity  profiles 
o  Option  for  mass  (vs  salinity)  flux 

Finer  control  over  energy  loan  ice  model 

o  Melting  point  can  be  linear  in  salinity 
o  Set  ice  minimum  and  maximum  thickness 
o  Set  ice  vertical  temperature  gradient 
o  Or  get  ice  surface  temperature  from 

New  diagnostics  within  HYCOM 

o  Time-averaged  fields  (in  archive  files) 
o  Synthetic  instrumentation 
o  3-D  particle  tracking 
o  surface  and  constant  depth  drifters 
o  isopycnic  drifters 
o  fixed  instruments  and  moorings 


HYCOM  2.2  (IV) 

Climatological  nesting  now  allowed 

o  Start  from  monthly  mean  outer  model  archive  files 
o  Allows  nested  runs  longer  than  the  outer  run 
o  But  with  less  accurate  boundary  state 
o  Probably  only  suitable  for  regional  nests 

Nesting  no  longer  requires  co-located  grids 

o  General  archive  to  archive  horizontal 
interpolation  (curvilinear) 

Hybrid  to  fixed  vertical  grid  remapper 

o  Allows  fixed-coordinate  nests  inside  hybrid 
coordinate  outer  domains 

o  HYCOM  to  (fixed-grid)  HYCOM 
o  HYCOM  to  NCOM 


HYCOM  2.2  (V) 

Enhanced  hycomproc  and  fieldproc 

o  NCAR-graphics  based 
o  Many  more  color  palette  options 
o  Can  read  in  an  arbitrary  palette 
o  Mark  locations,  and  draw  tracks,  on  plot 
o  Plot  diffusion  coefficients  and  tracers  (hycomproc) 
o  Overlay  vector  and  line-contours  (fieldproc) 

Added  fieldcell 

o  Like  fieldproc,  but  for  cell-array  (vs  contouring) 

o  Uses  NCAR’s  map  projections 

o  Typically  much  faster  than  fieldproc, 
but  can  leave  unfilled  cells 


ERA-15  SST:  Feb 


SSS  (hour  0} 


HYCOM  2.2  (VI) 

Diagnostic  fields  to  netCDF  and  other  file  formats 

o  Archive  fields  in  layer  space 
o  On  p-grid  (interpolated  velocity) 
o  3-D  archive  fields  interpolated  to  z-space 
o  On  p-grid,  or 

o  Sampled  along  arbitrary  tracks 
o  3-D  archive  fields  sampled  on  iso-therms 
o  Meridional  stream-function  from  (mean)  3-D  archive 
o  In  logical  array  space  (rectilinear  grids) 
o  Binned  to  latitude  bands  (curvilinear  grids) 
o  Atmospheric  forcing  input  fields 

o  Time  axis  depends  on  “.b”  file  format 
o  Any  “.a”  file  with  the  right  “.b”  file  structure 
can  be  converted  to  netCDF 

o  Fields  binned  into  lon-lat  cells 


HYCOM  CURVILINEAR  GRIDS  and  NetCDF 
Most  basin-scale  cases  use  a  Mercator  grid 

o  1-D  lat  &  Ion  axes  (rectilinear) 
o  Handled  well  by  many  netCDF  packages 

Global  HYCOM’s  Arctic  patch  grid  is  curvilinear 

HYCOM  netCDF  use  the  CF-1 .0  conventions,  which 
support  curvilinear  grids 

o  If  latitude  and  longitude  are  2-D  grids 
o  1  -D  axes  are  array  indexes 
o  Longitude  and  latitude  arrays  are  also  in  the 
file  and  identified  as  alternative  coordinates 

Most  netCDF  packages  are  not  CF-1 .0  aware 

o  Bin  into  uniform  lon-lat  cells  off-line 

o  Interpolate  to  a  1-D  latitude  and  longitude 
grid  off-line 

o  General  archive  to  archive  horizontal 
interpolation 

Archive  to  archive  remapper  can  also  be  used  for 
standard  (non-native)  grids 

o  MERSEA  grid  is  uniform  1/8° 
o  AOMIP  grid  is  rotated  uniform  1/2° 


GoM  NESTED  TEST  DOMAIN 
Same  resolution  nesting  unexpectedly  useful 

o  No  need  to  rerun  large  domain 
o  Change  atmospheric  forcing  (e.g.  use  MM5) 
o  Change  vertical  structure 
o  Tracer  studies  (e.g.  add  biology) 

1/12°:  Gulf  of  Mexico  inside  Atlantic 

o  Change  from  20m  to  5m  coastline 
o  Run  for  Aug  1999  to  equilibrate 
o  Run  Sep-Nov  as  standard  test  case 

Used  to  test  advection  schemes 

All  needed  file  are  prebuilt 

o  ftp://hycom.rsmas.miami.edii/awall/hycoin/GOMdO.08/ 
o  Uses  2.1 .20,  and  2.1 .34  is  also  available 
o  Includes  a  passive  tracer 

Also  tested  (new)  climatological  nesting 

o  Start  from  monthly  mean  outer  model  archive  files 
o  Climatological  nested  boundary  conditions 


DAILY  VS  MONTHY  NESTING 
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sea  surf,  height  Dec  01, 1999  OOZ  [06.1  H] 


CANDroATE  FEATURES  FOR  HYCOM  2.3 
Stable-code  vs  new  features 

o  Released  code-base  has  to  be  tested  and  stable 
o  New  features  can  be  a  significant  improvement 
o  Will  add  interim  releases  to  web  page 
o  Features  may  be  removed  in  next  released  code 

Fully  region-independent 

o  Compile  once,  run  on  any  region  and  any 
number  of  processors 

o  Run-time  memory  allocation 
o  Might  reduce  performance  (fewer  compiler 
optimizations  available) 

o  Needed  for  full  ESMF  compliance 
Improve  split-explicit  time  scheme 
Tidal  forcing 
Diurnal  heat  flux  cycle 
Equation  of  state  that  is  quadratic  in  salinity 
Wind  drag  coefficient  based  on  model  SST 
Initial  support  for  ESMF 


HYCOM  AND  ESMF 

Earth  System  Modeling  Framework 

http  ://w  w  w.esmf.ucar.edu/ 

o  Superstructure  couples  components 
o  Air/Ocean/Ice/Land 
o  Asynchronous  I/O  component 

•  Run  “concurent”  with  model  components 

o  Infrastructure  provides  data  structures 
and  utilities  for  building  scalable  models 

Add  a  superstructure  “cap”  to  HYCOM 

o  Simplifies  coupled  systems 
o  HYCOM  coupled  to  LANL  CICE  sea-ice 
o  Convert  atmospheric  field  processing  and  the 
energy-loan  ice  model  into  ESMF  components 

o  Use  ESMF  for  (user-level  asynchronous)  I/O 

o  Interoperate  with  other  ESMF  compliant  ocean 
models  (e.g.  Poseidon,  MITgcm,  MOM4) 

The  initial  ESMF  support  will  be  optional 

o  HYCOM  version  2.3 

ESMF  will  eventially  required  to  run  HYCOM 

o  HYCOM  version  3.0 


